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Abstract 

The high elastic modulus and interlayer strengths of graphite flakes make them a durable solid superlubricant. 
Apart from this, they have configurable electrical properties, exhibit quantum Hall effects, and possess a myriad of 
useful photonic properties. The self-retraction behavior of graphite flakes can have significant impact on the 
creation of ordered stacks for various applications because any accidental or intentional displacement of the top 
flake over the stacks below may result in a misalignment of the carbon-carbon atomic arrangement which, in turn, 
can have influence over the electrical and photonic properties. It has also been revealed that there was a tendency 
of the displaced microflake to fail at times to return to its original starting position and orientation. Here, we 
elucidate this behavior by considering the influence of the interlayer potential forces based on minimal potential 
energy seeking. The maps of the parameters interrogated here provide the ability for precautions to be 
undertaken. They also potentially permit the creation of an array of microflake stacks in which the metastable 
states permit different information to be encoded by virtue of the differentiated photonic or electrical 
characteristics readable from each array site. 



Introduction 

The many favorable mechanical, electrical, thermal, and 
biocompatible properties of graphite render it as an 
important material. Much work has been expended to 
investigate the superlubricity [1,2] between graphite 
layers where the high elastic moduli and interlayer 
strengths make it a durable solid lubricant. It also has 
configurable electrical properties [3-5], exhibits quantum 
Hall effects [6], and possesses a myriad of useful photo- 
nic properties [7-10]. Digressing to the mechanical 
property aspect of a closely related structure, the ability 
of nested shells in individual multi-walled carbon nano- 
tubes (MWCNTs) to slide has been known for some 
time [11,12]. These MWCNTs have comparable super- 
lubricity as graphite, in which the interwall shear 
strength against sliding ranges from 0.08 to 0.3 MPa. 
An interesting behavior of these extracted inner shells is 
in their ability to self-retract into the outer shells when 
the extracting force is removed [11]. This has been 
attributed to the action of van der Waals (vdW) interac- 
tion forces [11,12]. Such a feature has led to the 
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conception of MWCNT-based actuator bearings [13] as 
well as nanoelectromechanical oscillators which offer 
the possibility of operating at frequencies in the giga- 
hertz range [14,15]. More recently, a similar self-retrac- 
tion capability was demonstrated on two graphite 
microflakes [16] of equal sizes. Although not mentioned 
previously, such a behavior can have significant impact 
in the creation of ordered graphene stacks for various 
applications. This is because any inadvertent or inten- 
tional displacement of the top flake over the stacks 
below may result in a misalignment of the carbon-car- 
bon atomic arrangement which, in turn, can have influ- 
ence over the electrical and photonic properties. It has 
also been revealed that there was a tendency of the dis- 
placed microflake to fail at times to return to its original 
starting position and orientation [16] unlike the 
MWCNT. Here, we attempt to elucidate this behavior 
by considering the influence of the interlayer potential 
forces. 

Method 

The interlayer interaction between graphite layers has 
been investigated on perfect sheets [17-19] and on 
sheets with defects [20]. One method to model the 
interlayer interaction between graphite layers is to use a 
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Figure 1 The schematic description of the van der Waals interaction between the i-th and j-th atom. Schematic description of (a) the van 

der Waals interaction between the i-th and j-th atom in the registry-dependent potential approach. A depiction of (b) a top graphite flake (pink) 
that is moveable relative to a bottom (blue) graphite flake. 



Leonard-Jones type of potential However, such a model 
may not be sensitive enough to account for the registry 
of the honeycomb structures and its layers. It is gener- 
ally accepted that the interlayer interactions are domi- 
nated by long-ranged vdW interactions. A registry 
dependent interlayer interaction potential can be para- 
meterized (see Figure la) for layered carbon structures 
using [21,22] 

V(r ijA n iAnj ) = e-^-zo) [ C+ f{ Plj ) + f( Pji )] -a(^) " 

W m 

f^ = e -(p/8) 2 j:c2n(P/s) 2n w 

Pij 2 = r ij 2 - ( n i r ij) 2 

where V is registry dependent interlayer potential, 
subscript i indicates any atom on bottom flake, subscript 
j indicates any atom on top flake, r,y is distance between 
atom i and atom j, n is normal vector of corresponding 
flake, pij is transverse distance between two atoms, /is a 
function to estimate rapid decay of interlayer potential 
energy, z 0 = 3.34 A is equilibrium interlayer spacing, A 
= 3.629 A" 1 , C = 3.03, A = 10.238, S = 0.578, C 0 = 
15.71, C 2 = 12.29, and C 4 = 4.933. The inter-atomic dis- 
tance of carbon, a - 142 nm, provides a convenient 
depiction of geometrical size. In our simulations we 
apply a flake size of \\a and 6a in the x and y axis, 
respectively (see Figure lb). Figure 2 traces the potential 
variations as one perfect graphite layer over one of the 
same size as a function of the sliding distances in the 
xy-plane. It can be seen that there is a general increase 
in potential (more positive indicating less attraction) 
with sliding distance, thus, accounting for the tendency 
of the moving flake to self-retract when the motive 
force to slide is removed. However, the profile contains 
local potential peaks and troughs. Hence, while the 
motive force is strong enough to surmount these valleys 
in the process of moving the flakes apart, the retraction 



process may engender the flake ceasing to move at the 
first local trough encountered in the return trip. If we 
assume the mass of the flake to be small, it is possible 
to depict the retraction process as a series of quasi-static 
steps. In each of these steps, which should be small 
enough, the nearby landscape is interrogated in six ways 
Sx = Sy = ± 0.01 A, Sep = ± 0.5° to yield six potentials. 
In order to arrive at these values, we selected a few dis- 
placed positions and orientations and reduced the values 
until no changes in the return path. The flake can then 
be taken to advance to the incremental position or 
orientation based on the lowest of these potentials 
found if its value is lower than the current potential. 
Otherwise, the flake is assumed to have arrived at its 
final retracted state. 

In order to characterize the movement more systemati- 
cally, we adopt the following approach (see Figure 3 in 
which the distances apart are exaggerated for clarity). 
The top flake is originally at position I. It is then moved 
to position II some x-y distance away and rotated by an 
angle <fi to the #-axis. The flake is allowed to perform a 
series of potential searches at the vicinity of the center of 
the flake about ± Sx, ± Sy, and Sep to seek for the next 
position and orientation to retract to. This proceeds 
sequentially until the flake moves to a rest position III. 

Results and discussion 

The parameters tracked are (a) the center-to-center dis- 
tance between positions I and III (AC in Figure 3, (b) 
the angular displacement of the center-to-center posi- 
tions of I and III (a in Figure 4), (c) the center-to-center 
displacement between II and III (BC in Figure 5), and 
(d) the angular displacement at position III relative to 
the horizontal (0 in Figure 6). These parameters were 
traced in relation to x and y distance as well as ep intro- 
duced. Only one quadrant of possible movement was 
explored due to symmetry. 
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Figure 2 The potential distribution (units in meV) when a top graphite flake is moved. The potential distribution (units in meV) when a 
top graphite flake is moved relative to the bottom graphite flake of the same size (11c? versus 6a in the x and y axis, respectively, where a = 
0.142 nm) in the xy plane. The movement values are in terms of a. As movement should tend towards highest negative potential, the general 
landscape indicates a general tendency for the top flake to retract. However, the undulations imbue a possibility for it not to return to the exact 
original position and orientation. 




Figure 3 Depiction of the top flake orignally at position I that is moved to position II. Depiction of the top flake originally at position I 
that is moved to position II some x-y distance away and rotated by an angle 0 to the x-axis. The flake is then allowed to perform a series of 
potential searches at the vicinity of the center of the flake about ± Sx, ± Sy, and ± d<j> for sequential next position and orientation retraction that 
ends a rest position III. The parameters tracked are the center-to-center distance between positions I and III (AC), the angular displacement of 
the center-to-center positions of I and III {a), the center-to-center displacement between II and III (B,C), and the angular displacement at position 
III relative to the horizontal (0). 
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Figure 4 The map of the center-to-center distance between positions I and III. Map of the center-to-center distance between positions I 
and III (A,C) in Figure 3 at x-y distance away and rotated by an angle 4> to the x-axis. The inability to self-retract to the original starting position 
is evident. 



In the plot of distributions of AC using <fi = 0° in Fig- 
ure 4, we find that there is an ability for the top flake 
center to self-retract to a high degree back to its original 
position center (e.g. <a). There is a somewhat reducing 
ability to do this when x and y are increased. This is 
logical as there is a greater opportunity for the flake to 



encounter more local minima during retraction even 
though the overall potential difference is higher with 
greater displacement. Nevertheless, one finds an easier 
capability to self -retract fully when displaced in the y as 
opposed to the x direction. This demeanor is rather pre- 
served despite the angle <fi changing, as seen in the other 
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Figure 5 The map of the angular displacement of the center-to-center positions of I and III. Map of the angular displacement of the 
center-to-center positions of I and III {a) in Figure 3 at x-y distance away and rotated by an angle 4> to the x-axis. 
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Figure 6 The map of the center-to-center displacement between II and III. Map of the center-to-center displacement between II and III (B, 

C) in Figure 3 at x-y distance away and rotated by an angle 4> to the x-axis. 



plots with values ranging from 15° to 90° at intervals of 
15 . Overall, however, one finds that there are a myriad 
of states with which the top flake is able to stop at, 
depending on the displacing positions and orientations. 
The plot of distributions of a with <fi in Figure 5 essen- 
tially tells of the quadrant the flake returns to. There 
were cases wherein if AC was small, the flake may move 
into other quadrants. Nevertheless, we did not find 
unexpected cases in which the flake could move to 
another quadrant when AC was large. 

From the plot of distributions of BC with <j) = 0° in 
Figure 6, we find general rings which typically indicate 
that a flake that is pulled further apart will travel further 
to try to restore its original position. This is consistent 
with the potential map shown in Figure 2. Nevertheless 
one finds that there are islands when the flake is pulled 
more in the x-direction where BC is zero. This trend 
mirrors the finding given in Figure 4. From the plot of 
distributions of 0 with (p = 0° in Figure 7, we find that 
there will generally be large islands of non-rotated states 
(i.e. 9 9 0°). The other plots indicate that the top flake 
to mostly 9 to mostly rest at the angles with which 0 is 
used. This result infers that the retraction process over- 
all appears to seldom find S<p as the next lower potential 
step to advance to. However, the map also indicates that 
there are exceptions to this. We suspect that such cases 
should closely correspond when the top flake encoun- 
ters a peak in its path and, thus, undergoes rotation 
more in order to find the next potential minima. 



The landscape maps of Figures 4 to 7 provide us with 
the ability to analyze the trajectories of certain selected 
starting conditions. If we consider the case in which the 
top flake was moved initially to (x = 0.7a, y = 0.2a, (j) = 
0 ), we see from Figures 8 and 9 that it took slightly 
more than 50 steps before the flake became immobile 
over a short distance (about 0.2a). From Figure 10, the 
flake moved towards the x-axis (resulting in a = 0 ) 
when this occurred. Figure 11 informs us that the flake 
rotated when it encountered the local minima and then 
restored to its original orientation when it finally came 
to rest. The situation when the top flake was moved 
initially to (x = OAa, y = 0.2a, (j) = 0 ) is slightly differ- 
ent. While the flake encountered a local minimum to 
reside in, this happened at a longer distance away 
(about 3a) from Figures 8 and 9. The trajectories for a 
(Figure 10) and 9 (Figure 11) were nonetheless 
unchanged. A significantly larger behavior change 
occurred when the top flake was moved initially to (x - 
la, y - 0.2a, <fi = 30). With this, the movement trajec- 
tory for AC and BC were gradual before the flake even- 
tually came to a rest. The interesting behavior is 
denoted in Figure 11 wherein the flake is seen to rotate 
in order to overcome the potential barriers in the way 
of motion. This is the basis of lubricity as noted pre- 
viously [17]. When the top flake was moved initially to 
(x = 0.2a, y = 3a, (j) = 0 ), a significant degree of self- 
retraction was achieved (Figures 8 and 9). It is noted 
that the flake moved slightly towards the right (see 
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Figure 10) in the process but exhibited little rotation 
(Figure 11). 

We highlight certain caveats to the results provided 
here. The position and orientation tracing approach 
here assumes that there is no defect in either of the 
interacting flakes. The ability for graphite layers to be 
deform has been studied [23], and it is possible that 
such defects may strongly affect the self-retraction beha- 
vior [24]. This will be more evident for micron-scale 
flakes where the creation of entire interacting regions 



with perfect lattices is understandably more difficult. 
Such lattices should, however, be more achievable if the 
dimensions of the flakes are smaller. It is also important 
to note that the model here assumes that the force to 
move the top flake does not comprise an axial compo- 
nent that deforms the lattice. In previous demonstra- 
tions of self-retraction behavior, a mechanical probe 
that depressed and pulled the top flake laterally was uti- 
lized [16,24] in which such an effect cannot be ruled 
out. The measure of moving the flake with minimal 
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Figure 9 The trace of a. Trace of a in relation to step trajectories when the top flake was released at selected x-y distances away from the 
initial position and rotated by an angle 0 to the x-axis. 



axial deformation will most likely be accomplished using 
either optical [25], magnetic [26], electrical quadruple 
ions [27], or vibration [28] that have recently been 
reported with graphene sheets. The vibration of gra- 
phene nanostrips, it appears, offers the possibility of 
creating bending resonators with high sensitivity to 
environmental change [28]. It is conceivable, based on 
the findings here, for metastable misaligned states to 
appear with each cycle of actuation, to the extent of 
leaving a wrinkling effect on the sheets [29]. It is also 
noteworthy that long strips can result in gravity causing 
the structure to flop downwards. Remedy, however, is 
available by adopting a dangling arrangement, which has 



been demonstrated workable in applications associated 
with wetting monitoring [30,31]. 

Conclusions 

In summary, we have investigated the self-retraction of a 
graphite flake over another flake of the same size. This 
was done using a registry dependent interlayer interac- 
tion potential previously reported that is sensitive 
enough to account for the registry of the honeycomb 
structures and its layers. Under the assumption that the 
mass of the flake is small, the retraction process can be 
depicted as a series of quasi-static steps that seek the 
lowest potential. The results show that while there is an 
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Figure 10 The trace of BC. Trace of BC in relation to step trajectories when the top flake was released at selected x-y distances away from the 
initial position and rotated by an angle 4> to the x-axis. 



Ng et al. Nanoscale Research Letters 2012, 7:185 
http://www.nanoscalereslett.eom/content/7/1 /1 85 



Page 8 of 9 





7 0 -r 




60 - 




50 - 




40 - 










ft 


30 - 






■u 






20 - 




10 - 




0 - 




-10 - 



(x:7a ; y:G.2a ; t|3:3G 0 ) 



(x:7a ; y:0.2a : tt):0 a ) (x:4a : y:0.2a : c|):0 0 ) (x:0.2a ; y:3a^:0 0 ) 



o 



50 



100 150 200 2 50 300 3 50 400 450 500 

Step 

Figure 11 The trace of 9. Trace of 9 in relation to step trajectories when the top flake was released at selected x-y distances away from the 
initial position and rotated by an angle 4> to the x-axis. 



overall impetus to retract to the original position and 
orientation to restore to the lowest potential, there is a 
possibility for the return trajectory to encounter local 
potential minima that prevents the top flake from 
restoring fully. Essentially, this means that the graphite 
flake is able to assume meta-stable states. 

Such a behavior can have significant impact in the 
creation of ordered graphene stacks for various applica- 
tions. The maps of the parameters that we interrogated 
here provide the ability for precautions to be underta- 
ken. For instance, any movement in the y direction 
keeping x and <fi constrained will generally result in 
total retraction. Alternatively, it might be desirable to 
attain the meta-stable positions and orientations in 
order to obtain a different photonic or electrical beha- 
vior from an aligned stack. The maps then provide a 
road map to accomplish this. Such a behavior allows 
one to contemplate the creation of an array of micro- 
flake stacks in which the metastable states permit differ- 
ent information to be encoded by virtue of the 
differentiated photonic or electrical characteristics read- 
able from each array site. Using the maps, these meta- 
stable states may be altered to or 'erased' using the 
retraction behavior to the original stable state. This por- 
tends the possibility of high density data recording and 
retrieval at the nanometer scale. The understanding of 
the metastable states may also shed light into the wrink- 
ling behavior of graphene sheets, particularly if strips of 
them are to be used as resonators. 
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